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The equilibrium electric double layer (EDL) that surrounds the colloidal particles is determinant 
for the response of a suspension under a variety of static or alternating external fields. An ideal 
salt-free suspension is composed by the charged colloidal particles and the ionic countercharge 
released by the charging mechanism. The existing macroscopic theoretical models can be improved 
by incorporating different ionic effects usually neglected in previous mean-field approaches, which 
are based on the Poisson-Boltzmann equation (PB). The influence of the finite size of the ions seems 
to be quite promising because it has been shown to predict phenomena like charge reversal, which 
has been out of the scope of classical PB approximations. In this work we numerically obtain the 
surface electric potential and the counterions concentration profiles around a charged particle in a 
concentrated salt-free suspension corrected by the finite size of the counterions. The results show the 
large importance of such corrections for moderate to high particle charges at every particle volume 
fraction, specially, when a region of closest approach of the counterions to the particle surface 
is considered. We conclude that finite ion size considerations are obeyed for the development of 
new theoretical models to study nonequilibrium properties in concentrated colloidal suspensions, 
particularly the salt-free ones with small and highly charged particles. 



I. INTRODUCTION 

It is well-known that many transport phenomena in 
suspensions of charged particles, particularly their sta- 
bility and electrokinetic response, are closely related to 
the properties of the equilibrium electric double layer 
(EDL) surrounding the particles [1}(3]. Most of these the- 
oretical EDL models are based on the classical Poisson- 
Boltzmann equation (PB), a mean- field approach which 
has been the cornerstone for the study of the diffuse part 
of the EDL near charged interfaces in electrolyte solu- 
tions. 

Many of the studies concern suspensions with low par- 
ticle concentration, but nowadays it is the concentrated 
regime [1HZ| what deserves more attention because of 
its practical applications. From the theoretical point 
of view, these systems lack of a full understanding due 
to the inherent complexity associated with the increas- 
ing particle-particle electrohydrodynamic interactions as 
particle concentration grows. 

Very commonly, the electrical interaction between par- 
ticles is largely screened at typical electrolyte concen- 
trations, and particle ordering is impeded. However, 
if the concentration of external salt is sufhciently low, 
or there is not any salt added to the suspension, the 
electrical interactions can give rise to different short or 
long-ranged ordered phases at relatively low particle vol- 
ume fractions. These systems are usually called colloidal 
crystals or glasses [8l|9]. Suspensions composed just by 
the charged particles and their ionic countercharge (the 
so-called added counterions) in the liquid medium are 
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named salt-free suspensions, and its study has increas- 
ingly grown in recent years from experimental and theo- 
retical points of view. 

The theoretical models for such systems are generally 
scarce in comparison with those for suspensions in elec- 
trolyte solutions. In these systems there is a coupling be- 
tween the added counterions released by each particle and 
the opposite charge on the particle surface to preserve the 
electroneutrality. If the suspension is concentrated [S], 
the particle-particle electrohydrodynamic interactions in- 
troduce severe limitations when dealing with the response 
to external fields, not only for the mathematical difhcul- 
ties of dealing with many-body interactions, but also for 
the additional numerical problems raised when it comes 
to solve the equations, no matter the approximation cho- 
sen. The authors have addressed this problem by using a 
cell model approach to deal with electrokinetic and rhe- 
ological properties in concentrated suspensions [TUHIH]- 

In spite of the achievements of the cell model, it main- 
tains the drawbacks of the PB treatment, particularly, 
the neglecting of ionic correlations and the consideration 
of point-like ions in solution. Fortunately, for many sit- 
uations including low to moderate particle charges and 
monovalent ions in solution, this approach provides an 
excellent representation of the equilibrium problem. On 
the contrary, it has been argued that it yields unphys- 
ical ionic concentration profiles near highly charged in- 
terfaces, and also that it is not able to predict impor- 
tant phenomena like, for example, the overcharging, also 
called charge reversal or charge inversion 19 . In the lit- 
erature we can find different attempts to overcome these 
PB limitations. Some of them concerns microscopic de- 
scriptions of ion-ion correlations and the finite size of 
the ions [I9H23] . Others are based on macroscopic de- 
scriptions considering average interactions by mean-field 
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approximations, that include entropic contributions re- 
lated with the excluded volume effect when the ions have 
a finite size [^^HIO] . While the former ones are physically 
more robust to describe real systems, they are basically 
restricted to equilibrium conditions. The macroscopic 
approaches, however, permit us to make predictions un- 
der equilibrium and non-equilibrium conditions. Unfor- 
tunately, it seems that mean-field approaches eliminate 
some of the crucial information necessary for explaining, 
for example, the overcharging phenomena. In the past, it 
has been shown with Monte Carlo simulations, and also 
with some theoretical approximations like the HNC-MSA 
PT] , that electrical ion-ion correlations joined to exclude 
ion volume effects could be the key point to understand 
such phenomena [19) . 

In a series of papers, Lopez-Garcfa and coworkers [551 - 
[38] , studied the equilibrium EDL of a particle in a dilute 
colloidal suspension in an electrolyte solution account- 
ing for the finite size of the ions by using a modified 
Poisson-Boltzmann equation (MPB). They also assumed 
the existence of a common distance of closest approach 
to the particle surface for all ionic species. Their first re- 
sults indicated that the associated exclude volume effect 
was very important for equilibrium and non-equilibrium 
properties. Shortly after, their model was compared with 
some Monte Carlo simulations for general ionic condi- 
tions [HI, and while the latter predicted charge reversal 
for the case of multivalent ions, the MPB was unable to 
show such behaviours whatever the valence of the ions 
may be. It was concluded that the neglecting of ion-ion 
correlations in the MPB model was its main drawback. 

Quite recently, Lopez-Garcfa et al. [Hj have extended 
their previous model to allow for a different distance of 
closest approach to the particle surface for each ionic 
species. The most important fact is that it predicts 
charge reversal under appropriate circumstances by only 
considering such ionic excluded volume effects. To our 
knowledge, it is the first time that a phenomenological 
theory based on macroscopic descriptions for equilibrium 
and non-equilibrium properties in colloidal suspensions is 
able to predict charge inversion, mainly with multivalent 
ions. In other words, the charge reversal phenomenon, 
that has been elusive for such formalisms, now displays 
a reasonable explanation from the macroscopic approach 
based on mean-field approximations. No consideration of 
explicit ion-ion correlations is then required. Although 
no charge reversal is possible in an ideal salt-free suspen- 
sion because the whole ionic countercharge equals that 
of the particle, and cannot exceed it, these achievements 
support our choice of using a similar macroscopic ap- 
proach for our salt-free concentrated suspensions, with 
the help of a cell model. It has been recently shown that, 
up to moderately strong electrostatic couplings, the cell 
model proves accurate for the prediction of osmotic pres- 
sures of deionised suspensions in agreement with Monte 
Carlo simulations and renormalized-effective interaction 
approaches [i5] . 

The macroscopic approach provides reasonable predic- 



tions for many non-equilibrium phenomena in concen- 
trated suspensions, which are in good agreement with 
experiments. Other descriptions, like Monte Carlo's or 
microscopic models, are unable of making feasible pre- 
dictions out of equilibrium. Summarizing, in this work 
we aim at studying the EDL of a particle in a concen- 
trated salt-free suspension considering that the added 
counterions have finite ion size. Our approach is a mean- 
field one with an adequate MPB equation inside the cell. 
The procedure is a generalization of that already used by 
Borukhov (55] for the special case of a salt-free suspen- 
sion, valid for the concentrated case. Unlike Borukhov's 
treatment, our model also incorporates an excluded re- 
gion in contact with the particle of a hydrated radius size, 
which has been shown by Aranda-Rascon et al. [371 ISH] 
to yield a more realistic representation of the solid-liquid 
interface and also to predict results in better agreement 
with experimental electrokinetic data. The equations will 
be numerically solved and the equilibrium surface poten- 
tial and counterion concentration profiles inside a cell will 
be analyzed upon changing particle volume fraction, par- 
ticle surface charge density, and size of the counterions. 
In order to show the realm of the finite ion size effect, the 
results will be compared with MPB predictions that do 
not take into account a finite distance of closest approach, 
and also with standard PB predictions for point-like ions. 



II. THEORY 
A. The cell model 

To account for the interactions between particles in 
concentrated suspensions, a cell model is used (bare 
Coulomb interactions among particles are included in an 
average sense, but ions-induced interactions between par- 
ticles as well as ion-ion correlations, are ignored). For 
details about the cell model approach see the excellent 
review of Zholkovskij et al. 44J. Different electrokinetic 
and rheological phenomena in salt-free concentrated sus- 
pensions have been developed by the authors according 
to this concept: the DC electrokinetic response to static 
electric fields [TU] (electrophoretic mobility and electrical 
conductivity), their response to oscillating electric fields 
[TTl [12] (dynamic mobility, complex conductivity and di- 
electric response); the electroviscous effect TS^, and some 
related to more realistic descriptions of the phenomena 
taking into account the possibility of chemical reactions 
in the system in agreement with experimental conditions 
[imiSj . We have learned from those works how impor- 
tant the description of the EDL is for the non-equilibrium 
responses. 

Concerning the cell model. Fig. [T] each spherical par- 
ticle of radius a is surrounded by a concentric shell of the 
liquid medium, having an outer radius b such that the 
particle/cell volume ratio in the cell is equal to the par- 
ticle volume fraction throughout the entire suspension, 
that is 
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FIG. 1; Cell model including the distance of closest approach 
of the counterions to the particle surface. 



(2.1) 



The basic assumption of the cell model is that the 
macroscopic properties of a suspension can be obtained 
from appropriate averages of local properties in a unique 
cell. 



B. Finite size of the counterions 

Let us consider a spherical charged particle of radius 
a and surface charge density a immersed in a salt-free 
medium with only the presence of the added counterions 
of valence Zc- The axes of the spherical coordinate system 
(r, 9, ip) are fixed at the center of the particle. In the 
absence of any external field, the particle is surrounded 
by a spherically symmetrical charge distribution. 

Within a mean-field approximation, the total free en- 
ergy of the system, F = U — TS, can be written in terms 
of the equilibrium electric potential 5'(r) and the coun- 
terion concentration ricir). The configurational internal 
energy contribution U is 



U 



dr 



V*(r)p -I- Zcenc(r)^'(r) - ^lcncir) 



(2.2) 

The first term is the self-energy of the electric field, 
where eo is the vacuum permittivity, and is the rela- 
tive permittivity of the suspending medium. The next 
term is the electrostatic energy of the counterions in the 
electrostatic mean field, and the last term couples the 
system to a bulk reservoir, where fj,c is the chemical po- 
tential of the counterions. 

The entropic contribution —TS is 



TS = keTnl 
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where ks is Boltzmann's constant, T is the absolute tem- 
perature, and n™"^ is the maximum possible concentra- 
tion of counterions due to the excluded volume effect, 
defined as n™°^ = V~'^, where V is the average volume 
occupied by an ion in the solution. The first term inside 
the integral is the entropy of the counterions, and the 
second one is the entropy of the solvent molecules. This 
last term accounts for the ion size effect that modifies the 
classical Poisson-Boltzmann equation and was proposed 
earlier by Borukhov et al. [27] 

The variation of the free energy F = U — TS with 
respect to ^'(r) provides the Poisson equation 
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and the counterions concentration is obtained performing 
the variation of the free energy with respect to Ucir) 
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where be is an unknown coefficient that represents the 
ionic concentration where the electric potential is zero. 
Applying the spherical symmetry of the problem and 



combining eqns (2.4) and (2.5), we obtain the modified 



Poisson-Boltzmann equation (MPB) for the equilibrium 
electric potential 
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As the cell is electrically neutral we have 
Q — Aira^a — —An / p{r)r^dr 



AttZcB / (r)r dr 



(2.7) 



The unknown coefficient be can be calculated itera- 



tively substituting the value of nc{r), eqn (2.5), in the 
last expression 
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We need two boundary conditions to solve the MPB 
equation. The first one is 



dr 







(2.9) 
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which derives from the electroneutrality condition of the 
cell and the application of Gauss theorem to the outer 
surface of the cell. The second one is 



^/{b) = 



(2.10) 
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that fixes the origin of the electric potential at r = &. 



The MPB problem, eqns (2.6), (2.91 and (2.10), can be 



solved iteratively using eqn (2.8 1 to find the unknown co- 
efficient be- Fortunately, the authors have recently shown 
[l4j that this kind of problem can be solved avoiding the 
iterative procedure. For this purpose, we will use dimen- 
sionless variables, which are defined as 
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rewritting eqn (2.6 1 as 
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where we have defined the function g{x). If we differen- 
tiate it, after a little algebra, it is possible to eliminate 
the unknown coefficient and find that 

g'{x) + z,g{x)^'{x) + ^g\x)^'{x) = (2.13) 



where the prime stands for differentiation with respect 



to X. In terms of the electric potential, eqn (2.13) is 
rewritten as 



{x) + (x) 
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X 



= (2.14) 



Eqn ( |2.14 1 is a nonlinear third-order differential equa- 
tion that needs three boundary conditions to completely 
specify the solution. Two of them are provided by eqns 
(2.9) and (2.10), which now read 



*'(/i) = ^{h)=Q 



(2.15) 



where h — {b/a) ~ cj)^^/^ is the dimensionless outer ra- 
dius of the cell. The third one specifies the electrical state 
of the particle, and can be obtained by applying Gauss 
theorem to the outer side of the particle surface r = a 



d*(r) 



dr 

Its dimensionless form is 

*'(1) = -CT 



(T 
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Regarding dilute suspensions, it is very common in the 
literature to use the surface charge (also regulated sur- 
face charge) or the surface potential as a boundary condi- 
tion at the particle surface when solving the equilibrium 



Poisson-Boltzmann equation, and all of them are equally 
valid. When it comes to concentrated suspensions, we 
prefer the use of the particle surface charge as boundary 
condition by the following reasoning. The particle charge 
is a property that can be often measured experimentally. 
Instead, the surface potential depends on the choice of 
the potential origin, which is the solution "bulk" for the 
dilute case (see the discussion in the Appendix). As in 
the concentrated case there is not a clear "bulk" , the ori- 
gin is typically chosen at the outer surface of the cell, 
r = b, which depends on the particle volume fraction, 
eqn (2.1). In any case, we always can use the surface 
potential as a boundary condition, representing the po- 
tential difference between the particle surface and the 
outer surface of the cell. 

If we consider point-like counterions, n™"^^ = cx), eqn 
(2.14) becomes the expression obtained by Ruiz-Reina 
and Carrique \T^ . It can be easily demonstrated that 
solving the third-order problem is mathematically equiv- 
alent to finding the solution of the iterative problem, that 
is, that any function satisfying eqns (2.14), (2.151 and 



(2.17) also satisfies eqns (2.6| and (2.8) — (2.101, and vice 
versa. We have removed the unknown coefficient be us- 
ing this procedure and now the electric potential can be 
obtained numerically in one single step. 



Eqns (2.14), (2.151 and (2.17) form a boundary value 



problem that can be solved numerically using the MAT- 
LAB routine bvp4c |45] , that computes the solution with 
a finite difference method by the three-stage Lobatto IIIA 
formula. This is a collocation method that provides a 
solution that is fourth order uniformly accurate at all 
the mesh points. The resulting mesh is non-uniformly 
spaced out and has been chosen to fulfill the admitted 
error tolerance (taken as 10~^ for all the calculations). 

Once we have found the electric potenti al ^(x ) we can 
obtain the coefficient be by evaluating eqn (2.12) at x = h 



(2.18) 



The ionic concentration nc{r) is determined using eqns 



(2.5) and (2.11). 



C. Excluded region in contact with the particle 

Following the work of Aranda-Rascon et al. [37], we 
incorporate a distance of closest approach of the coun- 
terions to the particle surface, resulting from their finite 
size. We assume that counterions cannot come closer to 
the surface of the particle than their effective hydration 
radius, R, and, therefore, the ionic concentration will be 
zero in the region between the particle surface, r — a, 
and the spherical surface, r = a + R, defined by the 
counterion effective radius. This reasoning implies that 
counterions are considered as spheres of radius R with a 
point charge at its center. 
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The whole electric potential ^(r) is now determined 
by the stepwise equation 



r dr 



[Eqn (2.61 



a<r <a + R 
a + R<r <b 



(2.19) 



To completely specify the problem we now should force 
the potential to be continuous at the surface r — a + R, 
and also its first derivative, which is related to the con- 
tinuity of the normal component of the electric displace- 
ment at that surface, in addition to boundary conditions, 
eqns (2.9) and (2.10). Thus, in the region in contact with 
the particle [a, a + R], we are solving Laplace's equation, 
and, in the region [a + R, b], the MPB equation that we 



obtained previously in eqn (2.6) 



As we have seen before, we can eliminate the coefficient 
be changing the system of second order differential equa- 
tions, eqn (2.19), into one of third order, which hugely 



simplify the resolution process 



Eqn (2.14) 



1 < a: < 1 



1 + ^ < X <h 
(2.20) 

where we have use dimensionless variables. The bound- 
ary conditions needed to completely close the problem 
are 

' L\ ^ al (l + i?/a)2 

^p{h) = Q ^'p{h)=Q (2.21) 
^'i(l-|-f) = *p(l + f) ^^(l + f) = *'^(l + f) 

where subscript L refers to the region in which the poten- 
tial is calculated using Laplace's equation, and subscript 
P refers to the region in which we evaluate the MPB 
equation. The electric pote ntial i s obtained numerically 
in one single step using eqns ( 2.20 ) and ( 2.21 ). The be co- 
efficient and the ionic concentration 7ic(r) are determined 



by eqns (2.18) and (2.5), respectively. 



III. RESULTS AND DISCUSSION 

For all the calculations, the temperature T has been 
taken equal to 298.15 K and the relative electric permit- 
tivity of the suspending liquid = 78.55, which coincides 
with that of the deionized water, although no additional 
ions different to those stemming from the particles have 
been considered in the present model. Also, the valence 
of the added counterions has been chosen equal to 
+1 and the particle radius a — 100 nm. Other values 
for Zc could have been chosen. The model for point-like 
ions is able to work with any value of Zc in the Poisson- 
Boltzmann equation, but we think that the predictions of 
this model will be less accurate in the case of multivalent 



counterions, since it is based on a mean-field approach 
that does not consider ion-ion correlations. Nevertheless, 
when we take into account the finite size of the ions, the 
main objective of this work, we include correlations as- 
sociated with the ionic excluded volume, solving partly 
this problem, because we are still not considering the 
electrostatic ion-ion correlations. 

For the sake of simplicity, we assume that the average 
volume occupied by a counterion is ^ = (2-R)'^, being 2R 
the counterion effective diameter. With this considera- 
tion, the maximum possible concentration of counterions 
due to the excluded volume effect is n™"^ — {2R)~^. 
This corresponds to a cubic package (52% packing). In 
molar concentrations, the values used in the calculations, 
^max = 22, 4 and 1.7 M, correspond approximately to 
counterion effective diameters of 2R = 0.425, 0.75 and 1 
nm, respectively. These are typical hydrated ionic radii 
|46j . We think also that our model could be useful for 
short or medium length polyelectrolytes, using an effec- 
tive value for n™"^, as long as different ordered phases 
do not coexist. 

We will discuss the results obtained from three dif- 
ferent models, the classical Poisson-Boltzmann equation 
(PB), the modified Poisson-Boltzmann equation by the 
finite ion size effect (MPB), and the MPB equation in- 
cluding also the distance of closest approach of the coun- 
terions to the surface of the charged particle (MPBL) . 
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FIG. 2: Dimensionless electric potential distribution (a) and 
counterions concentration (b) along the cell for different ion 
sizes, considering (dashed lines) or not (solid lines) the ex- 
cluded region in contact with the particle. Black lines show 
the results of the PB equation. 
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Fig. [2] shows the dimensionless equihbrium electric po- 
tential distribution, Fig. [2k, and the counterions con- 
centration profiles. Fig. [2^^ along the cell. We display 
in solid black lines the predictions of the PB equation, 
coloured solid lines stand for the results of the MPB equa- 
tion and dashed lines correspond to the outcomes of the 
MPBL model. Different colours stand for different coun- 
terion sizes. The particle surface charge density have 
been chosen equal to —40 /xC/cm^ and the particle vol- 
ume fraction is (j) = 0.5 (very concentrated suspension), 
which implies a normalized cell size of b/a = 1.26. 

We can observe in Fig. [2^ that the inclusion of the 
finite ion size effect (MPB calculations) always rises the 
surface electric potential in comparison with the PB pre- 
dictions, which are recovered in the limit of point-like 
counterions (n™"^ = 00). The reason is that the limi- 
tation of the counterions concentration in the vicinity of 
the particle seriously diminishes the screening of the par- 
ticle charge, and consequently, there is an increment of 
the surface potential. An additional significative increase 
of the surface potential is obtained when we take into ac- 
count the distance of closest approach of the counterions 
to the particle surface (MPBL model). As expected, the 
existence of a Laplace region free of counterions also pe- 
nalizes the screening of the particle charge. 

We can also see in Fig. [2]d how the finite ion size 
(MPB) creates saturated regions which corresponds to 
plateaus in the counterionic concentration profiles. These 
condensates extend to larger distances when we consider 
the excluded region in contact with the particle (MPBL) . 

Fig. |3] shows the effective charge divided by the par- 
ticle charge along the cell for a particle volume fraction 
(/) = 0.5 and for different ion sizes, accounting or not for 
the excluded region in contact with the particle. The 
effective charge is the charge seen at a given radial dis- 
tance from the center of the particle which includes the 
own particle charge and the ionic volume charge density 
in the double layer integrated from the particle surface to 
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FIG. 3: Effective charge divided by the particle charge along 
the cell for different ion sizes, considering (dashed lines) or not 
(solid lines) the excluded region in contact with the particle. 
Black lines show the results of the PB equation. 
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FIG. 4: Dimensionless surface electric potential for different 
values of particle surface charge density (a) and particle vol- 
ume fraction (b) for different ion sizes, considering (dashed 
lines) or not (solid lines) the excluded region in contact with 
the particle. Black lines show the results of the PB equation. 



r distance [15]. From this definition, the ratio Qeff/Qp is 
equal to one at the surface of the particle and zero at the 
outer surface of the cell. The results of this figure rein- 
forced those discussed previously in Fig. [2j As expected, 
the larger the ion size, the softer the screening of the par- 
ticle charge. Besides, when we consider the distance of 
closest approach of the counterions to the particle sur- 
face, the screening is even less effective at every given 
distance from the latter. 

Fig. |4] displays the dimensionless equilibrium surface 
electric potential for a wide range of particle surface 
charge densities. Fig. |4^, and particle volume fractions. 
Fig. ^p. We repeat this study for different counterion 
sizes. We find that the surface electric potential increases 
with the particle charge density, Fig. |4^. However, in 
some cases there is a different behaviour in comparison 
with the point-like case. Initially, a fast and roughly in- 
crease of the surface potential with the surface charge 
density is observed, which is followed by a much slower 
growth at higher surface charge densities for the PB case 
or when the size of the counterions is very small. This 
phenomenon is related to the classical counterion conden- 
sation effect: for high surface charges a layer of counteri- 
ons develops very close to the particle surface [5]. When 
the ion size is taken into account (MPB) we limit the 
appearance of the classical condensation effect, because 
when the surface charge is increased, the additional coun- 
terions join the condensate enlarging it. If the region of 
closest approach of the counterions to the particle surface 
is also considered (MPBL), the mechanism is the same, 
but now the additional counterions are located in farther 
positions from the particle surface. This explains the fur- 
ther increase of the surface potential observed for large 
ion sizes and high surface charges densities. 

On the other hand, the surface electric potential de- 
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FIG. 5: Dimensionless surface potential against the surface 
charge density for different particle volume fraction values. 
Solid lines stand for the results of the PB equation. Dashed 
lines stand for the results of the MPBL model. 




FIG. 6: Dimensionless surface potential against the particle 
volume fraction for different surface charge density values. 
Solid lines stand for the results of the PB equation. Dashed 
lines stand for the results of the MPBL model. 



creases when the particle volume fraction increases, irre- 
spectively of the case studied: PB, MPB or MPBL, Fig. 
|4]d. When the particle concentration raises, the available 
space for the counterions inside the cell decreases and, 
consequently, the screening of the particle charge largely 
augments, thus reducing the value of the surface poten- 
tial. 

Figs. [5] and [6] expand the results of Fig. |4] In them, 
we compare the results of the PB equation (solid lines) 
with those of the MPBL model (dashed lines) for a given 
counterion size (typical of an hydronium ion in solution) . 
Fig. [5] presents the dimensionless surface electric poten- 
tial at a wide range of particle surface charge densities. 
The different coloured lines correspond to different par- 
ticle volume fractions. The most remarkable fact shown 
in Fig. [5] is the large influence of the finite ion size effect 



even for moderately low particle surface charge densities. 
While for PB predictions the surface potential hardly 
increases with surface charge for each volume fraction, 
the MPBL results display an outstanding growth for the 
same conditions, associated with the lower charge screen- 
ing ability of the counterions because of their finite size. 
This fact will surely have important consequences on the 
electrokinetic properties of such particles in concentrated 
salt-free suspensions, as has already been shown for dilute 
suspensions in electrolyte solutions by Aranda-Rascon et 
al. [371 

Finally, Fig. [6] shows the dimensionless electric po- 
tential at the particle surface against the particle volume 
fraction. The different coloured lines correspond to differ- 
ent negative particle surface charge densities. As previ- 
ously stated, the particle surface potential monotonously 
decreases with volume fraction at fixed particle charge 
density. For the counterion size chosen, the larger the 
surface charge, the larger the relative increase of the po- 
tential at every volume fraction. 

From the results, we think that it is clear that the in- 
fluence of the finite ion size effect on the EDL description 
cannot be neglected for many typical particle charges and 
volume fractions. 



IV. CONCLUSIONS 

In this work we have studied the influence of finite 
ion size corrections on the description of the equilibrium 
electric double layer of a spherical particle in a concen- 
trated salt-free suspension. The resulting model is based 
on a mean- field approach, due to its reasonable success 
in modeling electrokinetic and rheological properties of 
concentrated suspensions, particularly the salt-free ones 
which deserve a special theoretical interest. 

We have used a cell model approach to manage 
with particle-particle interactions, and modified Poisson- 
Boltzmann equations (MPB and MPBL) to account for 
such ion size effects. The theoretical procedure has fol- 
lowed that by Borukhov but with the additional in- 
clusion of a region of closest approach for counterions to 
the particle surface [33 [S^ . We also have made numer- 
ical calculations stressing and discussing the results for 
concentrated suspensions. The results have clearly shown 
that the finite ion size effect (MPB equation) is quite im- 
portant for moderate to high particle charges at every 
particle volume fraction, and even more if the distance 
of closest approach of the counterions to the particle sur- 
face is taken into account (MPBL model). 

This equilibrium model presented in the paper is on 
the base of future nonequilibrium models of the response 
of a suspension to external electric fields. Experimen- 
tal results concerning the DC electrophoretic mobility, 
dynamic electrophoretic mobility, electrical conductivity 
and dielectric response, should be compared with the pre- 
dictions of the latter models to test them. To perform 
such comparisons highly charged particles are required. 
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Existing highly charged sulfonated polystyrene latexes 
could be good candidates. 

As a final conclusion, we think that it is mandatory 
to apply these models to improve the existing theories 
for predicting non-equilibrium properties in concentrated 
suspensions. This task will be addressed by the authors 
in the near future. 




Appendix: Invariance of MPB equation under 
changes in the origin of the electric potential 

In this Appendix we demonstrate that the new MPB 
equation is invariant under changes in the origin of the 
electric potential, as it must be according to basic phys- 
ical grounds. This is not a trivial question, because the 
classical Poisson-Boltzmann equation for the low particle 
concentration case in electrolyte solutions is not invariant 
under changes in the origin of the electric potential: it is 
only valid when we take the potential zero at an infinite 
distance from the particle surface (bulk) . When we study 
salt-free and/or concentrated suspensions, we do not have 
a bulk, because the electroneutrality is achieved at the 
outer surface of the cell, and we have to take care about 
the correct invariance of the equation. In this case we do 
not have an "infinity" to set the origin of the potential 
and we have also an unknown coefficient be that must 
be correctly defined to account for the latter statements. 
Our PB description for the salt-free case pH] satisfies the 
above-mentioned invariance requirement. 

As we said before, the unknown coefficient be, that ap- 



pears in the MPB equation, eqn ( 2.6 ) , represents the ionic 



concentration where the equilibrium electric potential is 
zero. Our choice along this work has bee n ^ (6) — 0. Ac- 
cording to this definition and using eqn (2.5), we can see 
that nc{b) = be. 

The potential difference between two points r, ri does 
not depend on the origin of potentials. Fig. [7| 



*(r) - ^^(ri) = ^*{r) - 4'*(ri) 



(A-1) 



where the magnitudes without asterisks refer to the origin 
^'(6) — and those with asterisks, to an arbitrary origin, 
^'*(ri) = 0, for the electric potential. The relationship 
between the potentials referred to both choices of electric 
potential origins is 



ili(r) = **(r) -I- 'J'(ri) 



(A-2) 



Substituting the last expression in the MPB equation, 
eqn (2.6 1, we obtain 



d?-^*{r) ^ 2d4'*(r) 
dr^ r dr 



Zee 



be cxp ( ^— ) exp 



/_£ce*(ri; 



exp 



- 1 



FIG. 7: Changes in the origin of the electric potential. 



After a little algebra, and defining 
be 



I - fcTsfc? J exp 

it yields 

d2**(r) 2d**(r) 



2ce>I'(ri) 
fcflT 



(A-4) 



dr 



Zee 



^^exp(-^) 



Zee 



exp 



(r) 



eoCr 

where we have defined 



n*(r) 



(A-5) 



b* exD (-^^ 
"c exp ^ j,^. 



(r) 



1 



exp 



(r) 



- 1 



(A-6) 



Eqn ( A-5 ) shows that 5'*(r) satisfies the MPB equa- 
tion, eqn ( 2.6 ) , with a different coefficient 6* . Let's finally 



check that the concentration of counterions, n*{r), is the 
same irrespective of the origin of the electric potential: 
n*{r) — ne{r), in spite of the different values of the un- 
known coefficients be and 5*. 



By inserting eqn (A-4) in eqn (A-6), we obtain 

neir) 



(ri) 
T 



\ ( Zee* 

O^^PI — ^ 



exp 



(r) 



1 + 



be 



exp 



/ Zee** 



(r) 



- 1 



1 - 



exp 



Zce>I'(ri) 
fcflT 



be 



(r) 
T 



exp 



/fee* 



(r) 
T 



- 1 



ne{T) (A-7) 



(A-3) as we wanted to demonstrate. 
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